pde by direct integration method